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Abstract: Dynamics of collisionless plasma described by the Poisson- 
Vlasov equations is connected with the Hamiltonian motions of particles 
and their symmetries. The Poisson equation is obtained as a constraint 
arising from the gauge symmetries of particle dynamics. Variational 
derivative constrained by the Poisson equation is used to obtain reduced 
dynamical equations. Lie-Poisson reduction for the group of canonical 
diffeomorphisms gives the momentum- Vlasov equations. Plasma density 
is defined as the divergence of symplectic dual of momentum variables. 
This definition is also given a momentum map description. An alter- 
native formulation in momentum variables as a canonical Hamiltonian 
system with a quadratic Hamiltonian functional is described. A com- 
parison of one-dimensional plasma and two-dimensional incompressible 
fluid is presented. 
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1 Introduction 



The purpose of this series of papers is to study the geometric structures 
underlying both kinematical descriptions and dynamical formulations of 
plasma motion, and to provide a geometrical framework for the Poisson- 
Vlasov equations 

V 2 / (q) = -e J /(z) rf 3 p (1) 

f + £-V,/-eV^ / -V p / = (2) 

of plasma dynamics. These works are initiated from a detailed study of 
an unpublished review on plasma dynamics by Marsden and Ratiu [1] 
and inspired from the problems and ideas in Marsden and Morrison [2J. 
We refer to [3]- [5] for background materials and for references on early 
works such as [6J- [TD] - The Poisson-Vlasov system was first written in 
Hamiltonian form by Morrison [H]-[l6]. The Poisson bracket was then 
shown to be the Lie-Poisson bracket on the space of plasma densities 

HE-US!- 

The kinematical or Lagrangian description will provide us with the 
configuration space G = Dif f can (T*Q); the group of canonical diffeo- 
morphisms on the phase space T*Q of motions of individual plasma 
particles. We shall elaborate, in the next section, the geometry of 
Dif fcan{T*Q) which is the framework for the Hamiltonian (Lie-Poisson) 
structure of the Poisson-Vlasov equations. 

In section three, we describe the Poisson equation as a kinematical 
constraint on the dynamics of Eulerian variables. We show that the 
Poisson equation characterizes the set of zero values of the momentum 
map associated with the action of additive group of functions J-(Q) on 
the position space Q of particles. This is the gauge group of particle 
motion on the canonical phase space T*Q. 

Momentum map realization of the Poisson equation implies that the 
true configuration space for the Poisson-Vlasov dynamics must be the 
semi-direct product space J-"(Q) (§) Dif f can (T*Q) with the action of the 
additive group J~(Q) of functions given by fiber translation on T*Q and 
by composition on right with the canonical transformations. In order to 
adopt the configuration space Dif f can (T*Q), which has been customary 
in earlier treatments of the subject [5],[T7]-[T9], we rather proceed by 
adapting a constraint variational derivative. In doing so, we implicitly 
take the advantage of the facts that the Lie-Poisson structure on J r *{Q) 
is trivial and that the constraint is of first class. The dual vector space 
0* of the Lie algebra Q=3tham(T*Q) of Hamiltonian vector fields turns 
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out to be the space of non-closed one- form densities on T*Q. By sym- 
plectic duality, this can be identified with the space of non-Hamiltonian 
vector fields on T*Q. This gives the decomposition of the tangent space 
TT*Q = g©(g*)" on which one can start a geometric treatment for Leg- 
endre transformation [20] . One of the main theme of the present work is 
to introduce the dynamical equations on (g*)^— part of TT*Q and make 
their relations with the Vlasov-Poisson equations (CP and 02]) precise. 

In section four, we shall present the kinematical reduction of the 
dynamics on T*G. In the momentum coordinates (H, IP) of g* the 
Poisson-Vlasov equations take the form 

V 2 0n(q) = e J V g -n p (z)rf 3 p (3) 

^-* m (.,, + ^w 

<7TP (7) 1 

— A! = -x h (w ( Z )) - -5% (z) (5) 

at m 
and they admit Lie-Poisson Hamiltonian structure with a Hamiltonian 
function linear in momenta. 

Section five will be devoted to the usual density formulation of the 
Poisson-Vlasov equations. The reduced dynamics on g* has a further 
symmetry given by the action of the additive group J-(T* Q) of functions 
on T*Q. The momentum map g* — > F*(T*Q) = Den(T*Q) defines the 
plasma density function /. 

In section six, we remark that Eqs.fll]) and (jSJ) are not the only equa- 
tions in momentum variables leading to the Vlasov equation. We present 
another set of equations described by a canonical Hamiltonian structure 
with a quadratic Hamiltonian functional. 

In section seven, we compare one dimensional plasma with two- 
dimensional incompressible fluid. 

2 Motion of Collisionless Plasma 

2.1 Kinematical description 

We let QcK 3 denote the configuration space in which the plasma par- 
ticles move. The cotangent bundle T* Q is the corresponding momentum 
phase space. This has a natural symplectic structure given by the canon- 
ical two- form VLt*q = dq l A dpi, where we employ the summation over 
the repeated indices. In the sequel, will denote the natural iso- 

morphism T*T*Q — > TT*Q taking uniquely a one-form to a vector field 
on T*Q [2], [I]. We let t — > tp t be a curve such that for each t, <p t is a 
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canonical transformation of the particle phase space T* Q preserving the 
canonical symplectic two form Qt*q. Given a point Z £ T*Q regarded 
as a reference point or, a Lagrangian label, we let z = y^(Z), also writ- 
ten as z = (q, p) = (p(Z,t), denote the current phase space point or, 
the Eulerian coordinates of plasma particles. The phase space velocity 
is given by the time dependent vector 

z = ^(Z)=X,(^(Z))=X(z,t) (6) 

and it generates the flow ip t . Since ip t is canonical, X is infinitesimally 
Hamiltonian. We shall assume that it is globally Hamiltonian and write 
h(z,t) for the corresponding Hamiltonian function so that at each t, 
X = X h . 

2.2 Group of canonical diffeomorphisms 

The flow ip t of the Hamiltonian vector field Xh on T* Q is a one parameter 
family of elements of the group G = Dif f can (T*Q) of all transforma- 
tions of T*Q preserving the symplectic two-form Qj>*q. In this work, we 
restrict the discussion to the canonical transformations connected to the 
identity [3], [21]. <p t acts on left by evaluation on the space T*Q of refer- 
ence plasma configuration to produce the motion of particles. Thus, a 
configuration of plasma can be specified by an element of G. The right 
action of G commutes with the particle motions and constitute an infi- 
nite dimensional symmetry group of the kinematical description. This 
is the particle relabelling symmetry [22] . For the motion of particles on 
T*Q described by the left action 

G xT*Q — > T*Q : (<p, Z) -» L^(Z) = ip(Z) = z (7) 

the velocity field on T*Q is the vector X^ lying in the tangent space 

T V G = {X v : T*Q -> TT*Q\ r T * Q oX v = (p} (8) 

of G at if. Here, tt*q '■ TT*Q — > T*Q is the natural projection of the 
tangent bundle of particle phase space. That means, X^ is a map over 
the element if : T* Q — > T* Q of the configuration space. Vector fields on 
G are then defined to be the maps X : G — > TG whose value at <p 6 G is 
given by X v £ T^G. Since the motion is associated with the left action, 
the velocity field X is invariant under the right action : G — > G of G 
for all ip £ G. 

The space {X\ (R^.X = TR^ o lof 1 = X, £ G} of right 
invariant vector fields on G is isomorphic to the tangent space over the 
identity mapping of T* Q with the isomorphism given at each (p £ G by 

X h ^X h oip VX h eT td G (9) 
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where X h is the Hamiltonian vector field fl^ Q (dh) on T*Q whose flow 
is (p. Then, the Lie algebra 

= (X W (T*Q); -[, ]) (10) 

of G consists of Hamiltonian vector fields on T* Q and [ , ] denotes the 
standard Jacobi-Lie bracket, with conventions as in pE]. The isomor- 
phism in Eq.([9]) is just the relation 

X v {Z,t) = X h {z,t) = {X h o<p t )(Z) (11) 

between the Lagrangian and the Eulerian velocities at the point z. In 
coordinates, if we decompose a canonical transformation ip into position- 
momentum pairs as (p — (£, rj), then, we have 

av(z, t) = C ( z) A + 4, ( z) A = A + p,A = x ft( z, t) . (i2) 

By the identification X^ — > h modulo constants, and the homomor- 
phism 

[Xh,Xk] = —X{h,k} T * Q > (13) 

we have the identification g ~ (^(T*^); { , }t*q) of the Lie algebra 
with the space of functions on T* Q endowed with the canonical Poisson 
bracket. 

We define the dual vector space g* of the Lie algebra g=%ham(T* Q) 
to be the non-closed one-form densities on T*Q 

g* = {U id ® dfi G A\T*Q) ® Den(T*Q) : dn id ^ 0} (14) 

where cf/i G Den(T*Q) is a volume six-form on T*Q. At each point z G 
T*Q the space A 6 (T*Q) of six- forms is one- dimensional whose basis can 
be chosen to be the symplectic volume dfi = Qj>*q = ^t* q/\Qt* qAQt* q- 
With these definitions we have the non-degenerate pairing 

(X h , U id ®dfi)= < X h (z), U id (z) > dfi(x) 

Jt*q 

h(z) V, • nj d (z) d/x(z) 

T*Q 

= (h,v z -nj d ®dii) (is) 

of g and g* with respect to the L 2 -norm. Here, denotes the compo- 
nents of the vector T$ id = f^^IL^) = nj d (z) • and the pairing in the 
integrand is defined over the finite dimensional space T*Q. From the 
last line of the above equations we conclude that the pairing of algebra 
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with its dual is nondegenerate if V 2 • nj d (z) 7^ 0. Since, Qt*q is nonde- 
generate, this is equivalent to the condition dllid ^ 0. The definition of 
dual space g* of the Lie algebra g already implies that the identification 
of vector space g = 3tham{T* Q) with the space of functions IF{T*Q) can 
be extended to the identification of dual space with the space Den{T*Q) 
of densities on T*Q via the definition of density 

which does not vanish due to nondegeneracy restriciton in the definition 
of g*. That is, we have the one-sided correspondence 

U id ®dfiEQ* — Y f d/ie Den(T*Q) 

with the definition of / given by Eq. flTB]) (see also the internet supple- 
ment to J3]). We shall show later the precise way of obtaining density 
in the context of a momentum map. By the above identifications with 
function spaces the pairing between Lie algebra g=J r (T* Q) and its dual 
g* = Den(T*Q) takes the form of multiply-and-integrate 

{h,f®dfj)= [ h(z)f(z) d/i(z) . (17) 
Jt*q 

Remark 1 The configuration space Dif f can (T*Q) can also be given a 
description in terms of sections of the trivial bundle T*Q xT*Q — >T*Q 
where T*Qo is the particle phase space with Lagrangian coordinates Z 
and T*Q carries the Eulerian coordinates z. The total space T*Q xT*Q 
is symplectic with the two-form Q_ = £It*q — &t*q- A diffeomorphism 
<p '■ T*Q ^-T*Q is canonical if &t*q ~ ^>*VLt*q — 0. It follows that 
Q_ vanishes when restricted to the graphs of canonical diffeomorphisms. 
Graphs are elements of the space T(T*QqxT*Q) of sections of the trivial 
bundle T* Q xT* Q -^T*Q . For a base point Z eT*Q , the total space is 
twelve dimensional and the graph (Z, cp (Z)) of a diffeomorphism is a six 
dimensional subspace. When if is canonical, f2_ vanishes on graphs and 
such graphs are called Lagrangian subspaces. If we denote the space of 
all sections of the trivial bundle on which the restriction ofQ- vanishes, 
namely the space of all Lagrangian sections, by Lag(T(T* Q xT* Q)) , 
then we have the identification 

Diffcan(T*Q) ~ Lag(T(T*Q xT*Q)) . 

To have the corresponding description for the Lie algebra of Hamiltonian 
vector fields, first observe that we can identify the dual space q* with the 
space of vector fields 

(g*) s = { ILL = 4*q(IW) e TT*Q | dli ld £ } 
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which are not Hamiltonian by definition. Then, we have the decompo- 
sition of the tangent space TT*Q = g©(g*)" into the underlying vector 
spaces of the Lie algebra and its dual. It turns out that the Lie alge- 
bra of Hamiltonian vector fields can be identified with the space of all 
Lagrangian submanifolds ofT(TT*Q) with respect to the Tulczyjew sym- 
plectic structure on TT*Q f2&^ . 

To extend the pairing of g and g* to a pairing of tangent and cotan- 
gent bundles we define an element 11^ of the covector space A^G of 
one-forms as a map over ip 

A\G = {n^ : T*Q ^ T*T*Q , vr T * Q oU^ = p} (18) 

where ht*q '■ T*T*Q -^T*Q is the natural projection, and a 1-form field 
II on G to be a section II : G — > A 1 G such that for each p G G we have 
G A^G. The cotangent space at p is then the space 

T;G={Il v ®d}xeAlG®Den(T*Q)} (19) 

of one-form densities on G = Dif f can {T* Q). The pairing of tangent and 
cotangent spaces at p G G becomes 

(x v , n v ®d/i)= [ < x v (z), iyz) > <^(z) . (20) 

Jt*q 

Following diagram summarizes the mapping properties with reference to 
particle phase space of elements of T V G, T*G, g and g* 



I I Q f^ Q , T*T*Q 

Xip ^ x h tir T * Q 7r T , Q ;tn id ^ Uip ■ 

T*Q p } T*Q = T*Q ( <p T*Q 

The adjoint action Ad v : g — > g of G on its Lie algebra is given by 
push forward of Hamiltonian vector fields on T* Q 

Ad v (X h ) = Tp o X h o p- 1 = p^(X h ) . (21) 

The coadjoint action (Ad v )* : g* ->■ g* of G on g* is the dual Ad* x with 
respect to the pairing in Eq. ffTo"]) of the map Ad^-i and is given by push 
forward as well 

Ad*! (n«) = ^ (JUd) = T* d (R v o L v -i) o U ld . (22) 
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Taking the derivatives at the identity one finds that the adjoint action 
ad Xh of g on itself and the coadjoint action ad* x of g on g* are generated 
by the Lie derivative with respect to the Hamiltonian vector field Xh G T 
T*Q. In other words, given Xh,Xk 6 g and IL^ e g*, the Lie derivative 
Cx h may be regarded as tangent vectors adx h £ Tx k g or adV € Tu id Q* 
depending on whether it is associated with adjoint or coadjoint actions, 
respectively 

The cotangent lift to T^G of the right action is T* or ^_ 1 R^(Il v ) = 

Hp o ip~ . For tp = if this gives the translation of the one form LT^ at if 
to the identity which we denote by Hid = Hp o f~ l . By definition, this 
is an element of the dual space g* of the Lie algebra g of Hamiltonian 
vector fields. The right invariant momentum map Ji : T*G — > g* for the 
lifted left action (i.e. the plasma motion) is defined by 

(Sl (n v ) <8> dfi, X h ) = (Hp ® d/i, T id Rp o X h ) 

= (n«<8)^,X ft ) (23) 

so that we have Ji^II^) = Il id G g* |]. 

3 Poisson Equation as a Constraint 

3.1 Momentum map description of Poisson equa- 
tion 

The canonical symplectic structure on T* Q is invariant under translation 
of fiber variable by an exact one-form over Q. This is the gauge trans- 
formation of canonical Hamiltonian formalism. If we identify T*Q with 
the space of one-forms A X (Q) on Q, this invariance may be described as 
the Hamiltonian action on T*Q 

A°(Q) x A\Q) -> A\Q) : (0 (q) , p-dq) -> p dq + dcf> (q) (24) 

of the space A°(<2) = J-(Q) of zero- forms. From an algebraic point of 
view, the exterior derivative d : A°(Q) — >■ A X (Q) can be interpreted as 
a map describing a Lie algebra isomorphism of the additive algebra of 
functions J~(Q) into the additive algebra of one-forms A l (Q) (23]. As 
the dual of any Lie algebra isomorphism is a momentum map, we have 

Jf(q) : A 2 (Q) Den(Q) (25) 

where the space of two forms A 2 (<2) on Q and the space of densities 
(three- forms) on Q are the duals, with respect to the L 2 — norm, of A 1 (Q) 
and A°(<2), respectively. We can identify A X (Q) and its dual space A 2 (Q) 
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by the Hodge duality oparator * associated with a Riemannian metric 
on Q. Then 

(Xf(Q) (p-dq, *# (q)) , (q)> = - / (q) d * d0 (q) (26) 

where *g?0 (q) is considered as a two-form over the one-form pdq, and 
we used the pairing of A 1 (Q) and A 2 (Q) given by integration over Q. 
Thus, the momentum map is 

Xf(Q) (p-dq, *d(j) (q)) = —d * d<j) (q) G Den(Q) . 

For the Euclidean metric on Q, the operator g? * d is the usual Laplacian 
Vg in Cartesian coordinates. 

The action of the additive group F(Q) of functions can be carried 
over objects defined on T*Q such as the group of canonical diffeomor- 
phisms, its Lie algebra and the dual of the Lie algebra. For our purpose of 
obtaining the Poisson equation as a constraint described by a momentum 
map, we restrict ourselves to the action of J-{Q) on the space J : {T*Q) of 
functions on T*Q. It will be convenient to describe the momentum map 
as the dual of some Lie algebra isomorphism into. We think of J-(T*Q) 
equipped with the Poisson bracket to be an algebra isomorphic to the 
Lie algebra of Hamiltonian vector fields (c.f. section [2721) . Then, F(Q) 
is a commutative subalgebra of (J-"(T*<2), { , }t*q) corresponding to the 
generators of the action on T*Q by fiber translation p i— )■ p — V q <p (q). 
The Hamiltonian function is — <p (q) £ J~"(Q.) C J-(T*Q). Thus, the re- 
quired Lie algebra isomorphism is from the additive algebra of functions 
J~(Q) into the Poisson bracket algebra on J-(T*Q). Together with the 
dualization we have 

(HQ),+) ► (HT*Q),{, h* Q ) 

t t 

Den(Q) < J ^ (g) Den(T*Q) 

and the momentum map Xf(q) : Den(T*Q) — > Den(Q) is computed 
from 

<Jf(Q) (/ (z) d/i (z)) , (q)> = - / / (z) (q) dfi (z) (27) 

Jt*q 

to be the volume density 

I HQ) (f (z) dp (z)) = -( / / (z) rf 3 p) d 3 q (28) 
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on Q [TB]. Combining the actions of J'(Q) on A 1 (Q) and F(T*Q) we 
have the momentum map 

J-F(Q) : A 2 (Q) x Den(T*Q) — > Den(Q) 

given by 

J-F(Q) (*# (q) , e/ (z) rf/i (z) ; (q)) = -(V 2 (q) + e J / (z) d 3 p) d\ 

whose zero value is the Poisson equation. This constraints the region in 
the product space A 2 (Q) x Den(T*Q) for consideration of dynamics in 
the Eulerian variables (0, /), namely 

^( Q )(°)/^(Q) = {(*d(j)Jdn) G A 2 (Q) x Den(T*Q) | 

V 2 0(q) + ej7(z) d 3 p=0} 

is the reduced space for the Lie-Poisson description. This corresponds to 
a subset of the dual space A 2 (<2) x g* of the Lie algebra A 1 (Q) x g with 
trivial bracket on the first factor. The underlying Lie group is J~(Q) x 
Diffcan(T*Q) with the first factor acting on canonical diffeomorphisms 
by composition with fiber translations. 

3.2 Constraint variational derivative 

The momentum map description of Poisson equation implies that we 
have to consider the configuration variables of collisionless plasma mo- 
tion to be (<f), (p) where the electrostatic potential is a function on 
Q and (p is a canonical diffeomorphism of T*Q generating the parti- 
cle motion. Hence, the configuration space of plasma motion must be 
J r (Q)®Dif f can {T* Q) where (§) denotes the semidirect product of groups 
with the additive group J~(Q) of functions acting on the second factor by 
composition on right. Here, we want to adopt an approach allowing the 
possibility to use much simpler configuration space G = Diff can (T*Q). 
More precisely, we want to use the Poisson equation as a constraint for 
the variational derivatives of Eulerian variables, in particular, the plasma 
density function. 

To this end, we consider the Green's function solution 

<f> f (q,t) = e[ K(q|q)/(z) d»(z) (29) 
Jt*q 

of the Poisson equation ([1]) which relates the plasma density / and the 
electrostatic potential (pj [H] , [H] , [25] . As an example of Eulerian quan- 
tities, we take the Hamiltonian function 
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H LP (f)= [ f(z)h f (z)dfi(z) (30) 
Jt-q 

of the Lie-Poisson formulation [H] , [E] , [H] , [25] . Here, the density de- 
pendent function 

hf (*) = 7^ + (31) 

is related to the Hamiltonian function h governing the particle dynam- 
ics up to a multiplicative factor in potential term. The Hamiltonian 
functional in Eq. (l30l is the total energy of the plasma in Eulerian coor- 
dinates. 

Lemma 2 Jiffi For the functional in Eq. / fffOj) with Eq. I[3l\) we have 

-p 2 + e0 / (q) = h(z). 



SH LP (f) _ 1 2 



5f 2m 1 

Proof. Using the Poisson equation and the Green's function solution, 
Hipif) can be put into the form 



H LP (f) = J^ ^ P 2 f(zW(z) 

□ 2 



2 



+ C - I I /(z)K(q|q)/(z) d,x(z) dfi(A) (32) 

IT*Q JT*Q 



up to the integral of the divergence term V q ■ (<^(q)V g 0y(q)). It is now 
easy to obtain the lemma where a factor of 2 comes from the symmetry 
of the Green's function f77] / ; f7^]/, fg^/ . ■ 

The constraint imposed by the Poisson equation is, in the language 
of Dirac formalism, first class and hence does not affect the Poisson 
bracket on the reduced space [26]. Thus, in obtaining equivalent dy- 
namical formulations in alternative Eulerian variables we must use the 
same constraint. The foremost example of such a variable arises from the 
identification of the dual space q* of the algebra of Hamiltonian vector 
fields with the space of densities Den(T* Q). In the more basic formula- 
tion of dynamics with the momentum variables G q* the Hamiltonian 
functional turns out to be 

H LP (U id ) = [ (Tl ld (z) , X hf (z)) dfi (z) (33) 
Jt*q 

which is equivalent to the functional H.Lp{f) under the identification 

(TH. 



Lemma 3 For the Hamiltonian functional in Eq. (E3|) we have 

-lfk^ = X " {Zl 
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4 Momentum Formulation of Dynamics 



4.1 Lie-Poisson dynamics 

We apply the standard Lie-Poisson reduction to G [Tf] , [18] , [27] . The 
right invariant extensions to T*G of elements of q* are obtained through 
composition with the momentum map 1^(11^) = IL^. In particular, for 
a functional H : g* — > K. the right invariant extension is the functional 
H R : T*G — > R defined by 

h r = h o j L , ^ n v ) = #(n id ) = h (n^ o c^ 1 ) 

and applying the chain rule with 11^ = Ui d o if we have the differential 
r ttr 5H r , ^SH r SH SU id 5H su id 

If n id (z) = U a (z)dz a = rL»dg' + IT(z)^ we have (<ffl id /5iy \ id = 1 
and 

|m= -^-U a (z)dz a \ id = -dll id (z) 

with ITjd denoting the components of the one-form 11^. Thus, for the 
differentiation of the right-invariant functionals we find 

SH R 5H 5H R 5H 

m; UG -m- d ' -^ lid ° = - { M- d ' Vz)Uid (34) 

where the overhead arrow denotes the components of Lie algebra ele- 
ment. One can now evaluate the canonical Poisson bracket on T*G at 
the identity using the above relations. This gives (+)Lie-Poisson bracket 
on g*, that is, {F R ,G R } T . G ^ {F, G} +LP . 

Proposition 4 Let U id G g* and [ , ] be the Jacobi-Lie bracket on g. 
Then the Lie-Poisson bracket on g* is given by 



{H(U id ),K(U ld )} LP = [ U id (z 



5H 5K 



m ld (z) 



dfx(z) (35) 



where 5H/5Ui d is regarded to be an element of g. 



The derivation can also be found in the internet supplement to [3]. 
The Lie-Poisson dynamics on g* may be written as follows: 



12 



Proposition 5 The Hamiltonian vector fields on g* for the Lie-Poisson 
structure defined by the bracket in Eq. < fff5]j have the form 

= -ad* SH/su . d (U u ) = C SH / S n id (Ili d ) = J LP (Ili d )-^- (36) 

where the Hamiltonian operator defining the bracket (ESP is given by 

JLp(Rid) = I n .J_ , % . jjj YlijL + jt.]jj) ( 37 ) 

\ l dp j dq 1 dpj dpi J 

w%th ~rj • n ? - = ~7T" ~\~ n ? ~7 etc. 

dq l J dq J dq l 



The operator in Eq. fl37|) is appearently skew adjoint with respect to 
the L 2 -norm and satisfies the Jacobi identity by construction 
The relation of JLp(Jl id ) to the Lie-Poisson bracket in Eq.( 13"5|) is 



n 



id 



5H 5K 



SU id 5U id 



5H , . . 8K . „ 8K . 5H , 

^)(^)-v,-— (n id .-—) 



<$n id 5n id 5n id <5n id 

where the divergence term on the right disappears upon integration. 

Remark 6 The Hamiltonian operator JipiJ^-id) may be considered to 
be a map taking a Lie algebra element Xh in g to the corresponding 
generator ad* Xh of the coadjoint action of g on q* . That is, 

Ji,p(JUd) : g — ► (-ad* : g* -»• g*) 

and with reference to particle phase space we have 

J LP : TT*Q — ► TO C Tg* C TT*T*Q 

where 0(Jli d ) = (Ad^-i(U.id) = (p t (Hid) | if £ C?| zs i/ie coadjoint orbit 
through Ui d G 0* . 

4.2 Vlasov equations in momentum variables 

In the Lie-Poisson setting, the Poisson- Vlasov equations arise from the 
Hamiltonian functional that generates the kinematical symmetries. This 
is the momentum function defined by means of the momentum map 

J L (X k ) (U id ) = [ (U id (z) , X k (z)> dfi (z) (38) 
Jt*q 

for the lifted left action of G. Here the generator X k G g is yet to be 
specified. Due to the constraint imposed by the Poisson equation, that 
is the constraint variational derivative, the relevant Hamiltonian vector 
field X k in Eq. (J3"8"]) turns out to be the one associated with the function 
hf given in Eq. ()3"T]) . 
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Proposition 7 For the right invariant Hamiltonian functional in Eg. / fffgj) 
the Lie-Poisson equations on g* are 



dUj (z) 



V h (IP(z)) + e 



gV/ (q) 

dqidqi 



(39) 



dIP (z) 

eft 



^(IT(z)) 



m 



1 




(40) 



with the constraint 



V^ n (q) = e / V 9 ■ n p (z) d 3 



p 



(41) 



Eqs. (I39j) and (|40p . which we call the momentum- Vlasov equations, 
follow from the Lie-Poisson structure expressed in coordinates of the mo- 
mentum one-form IT^ and will be shown, in the next section, to give rise 
to the Vlasov equation in the density variable. Phe proof of proposition 
follows from the constraint variational derivative of HLp(J\ i( i) obtained 
before and from Eq. (l36|) by computing, for example, the Lie derivative 
of Hid with respect to X h . 

5 Density Formulation of Dynamics 

Phe definition of plasma density is motivated by the following observa- 
tion. Regarding the variational derivative of functions on g* as elements 
of g means that there are functions h, k on T* Q such that the Jacobi-Lie 
bracket in Eq. (l35|) can be written as canonical Poisson bracket of h,k. 
Phen the Lie-Poisson bracket on g* becomes 



which requires, as in the definition of dual algebra, the divergence of the 
vector Ill d to be non-zero. 

5.1 Introducing the plasma density function 

We first show that the definition in Eq. (|T6l) of plasma density function 
leads to the correct Lie-Poisson structure in this variable. Phen we 
show that the relation between formulations of plasma dynamics in the 
momentum variables PL^ and the plasma density function / can be made 
precise in terms of a momentum map. We shall remark, in the next 
section, by presenting a canonical Hamiltonian system in momentum 
variables, that Eqs. (!39l and fT40l are not unique in the sense that they 
yield the Vlasov equation in density variable. 



<9ip(z) dlP(z) 

dpi dq l 



) {h(z),k(z)} T * Q dii(z) 



(42) 
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Proposition 8 The Hamiltonian operator Jlp (Hid) on g* transforms 
into the Hamiltonian operator 

JlpU) = V p / • V, - V g f ■ V p (43) 
on the space Den(T*Q) of densities under the correspondence in Eq. [To]) . 

Proof. Regarding the definition f lT6|) as a transformation of Eulerian 
variables we can compute the tranformation of the Hamiltonian operator 
J lp (Iljd) as follows. The derivative of / in the direction of Hid is a 3 x 6 
matrix of differential operators 

D/(IWz)) = [(A) -(A)] = [V P -VJ (44) 

which transforms the Hamiltonian operator Ju^n^) according to 

J LP (f) = D f Ol id ) ■ J LP (Tl id ) ■ D}(H id ) (45) 

where is the adjoint of Df with respect to the L 2 — norm [29J. A 
direct computation from Eq.( |4"5j) gives the operator Jz,p(/). ■ 

Remark 9 Dj transforms the Lie algebra elements, that is if H is a 
functional on Den(T*Q), then 5H/5f is a function onT*Q andD*A8H/5f) 
gives the components of the Hamiltonian vector field in q corresponding 
to the Hamiltonian function SH/5f. 

Remark 10 It can also be verified that the momentum- Vlasov equations 
and fiJW yield the Vlasov equation (0|) for f from the definition of 
the density. 

We recall that the Lie-Poisson structure of the Vlasov equation in 
density variable is defined by the bracket 

{H(f),K(f)} LP = [ /(z) \^- v ^t) dfi(z) (46) 
Jt*q 1<V(z) <V( z )Jt*q 

associated with the operator JlpU) and the Hamiltonian function H L p(f) 

5.2 Momentum map description of density 

Proposition 11 The momentum map for the action of the additive 
group J-(T*Q) of functions on g* defines the plasma density. 
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Proof. By definition in Eq. fll4j) of g*, the one-form IT^ is invariant 
under the addition of an exact one-form on T*Q. So, J r (T*Q) is the 
gauge group in the definition of q*. With reference to the identification 
of q with J r (T*Q), the action of J r (T*Q) on q* by translation can be 
regarded as the action of g on q* by 

(X k , U id ) ^ U ld + Vt T , Q (X k ) (47) 

where = (f^„ Q ) _1 . This can be interpreted as the action of the 

underlying vector space of g. Thus, we consider the Lie algebra isomor- 
phism F(T*Q) — > Q:k — >X k for the gauge equivalent classes of one-forms 
in q*. The dual of this is the required momentum map 

Jtr : 0* — > F*(T*Q) = Den(T*Q) (48) 

for the definition of the plasma density from the momentum variables. 
From definitions, we have 

(J tr (n id ) , k) = (U id , X k ) = (v z oQ^ Q oU id ,k^ (49) 

where V z is taken to be the dual of the exterior derivative d. When eval- 
uated in Eulerian coordinates z the momentum map ( 1491) gives exactly 
the definition (fi~6l) of the plasma density function. ■ 

Remark 12 It has been argued that the physical initial conditions must 
satisfy /(z,0) > [Lfflj fSUj . This requires the description of density 
by elements 11^ G g* with V z ■ nj d (z) > 0. Equivalently, in the lan- 
guage of differential forms, we have d(Jli d A Q^ tQ ) > 0. Consider a six 
dimensional domain D in T*Q with boundary dD. Then, the positive 
divergence implies 

[ Il id {z) A fl|. c (z) > (50) 

JdD 

so that we have a volume element or an orientation for the five dimen- 
sional boundary of the region D . This can now be related to the non- 
degeneracy of the orbit symplectic structure on q* . An element of the 
orbit through Ui d will be of the form Cx k (n^) ■ By definition, the orbit 
symplectic structure is 

iln id (C Xk (U id ),ji: Xg (U id ))= [ {g(z),k(z)}U ld (z)A Q 2 T , Q (z) 

JdD 

which, by Eg. ( T5flj) does not vanish for arbitrary functions g and k. 
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6 Equivalence of Momentum and Density Formu- 
lations 

Proposition 13 H LP {f) = H LP (U id ) 

Proof. Replace / in Eq.( l32|) by its definition to get 

HlpU) = I t^p 2 (v p ■ n,(z) - v, ■ n p (z)) d»{z) 

Jt*q ^ m 



e 2 
+ 2 



J J V q ■ n p (z)K(q|q)V 9 > ■ n p ,(z) rf/i(z)rf/i(z) 



T*Q 

upto divergence. The first and the second integrals are equivalent to 

2 

- P • n,(z) , - i-n p (z) • V,( K(q\q)V g/ ■ U^i) ), (51) 
m 2 

respectively. Using Green's function solution we obtain H LP (Jl id ). Con- 
versely, starting from the function H LP (IT^) we compute 

Hlp {Tkd) = [ (-V p /i/(z) ■ n,(z) + V q h f {z) ■ U p (z)) dfi (z) 
h f( z )f( z ) d ^( z ) 

T*Q 

which verifies the equivalence of the Hamiltonian functionals of the Lie- 
Poisson structures in momentum and density formulations ■ 

Proposition 14 The canonical Hamiltonian system with respect to the 
symplectic two form 



(U i ,U i )= [ (z) A <&n* (z) d/x(z) (52) 

Jt*q 

and for the Hamiltonian functional 

= J T , Q ( n ^< n '> + L s " u ' u ' + iH? 1 ™) (z) * (z) 

(53) 

which is quadratic in the fiber coordinates of T*T*Q gives the Vlasov 
equation. 

Proof. First of all, the density "Ho( z ) of the Hamiltonian functional Hq 
satisfies the divergence equation 

<)H,,{Z) - V z ■ ({Ui(X h (&) + -a%))X^ (z) = (54) 



at \ m 
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which is the conservation law in Eulerian form. To obtain the canonical 
equations for the Hamiltonian functional in Eg. [blfy we substitute the 
expressions for Xh and<ftj into Hamiltonian functional Hq and rearrange 
the terms to obtain 

h (u jt if) = / n, (z) + ^n,) (z) dp (z) 

2 ff , . , dW (z) d 2 K (q, q) 

-e yyn j( z)n(z)^- ^— {*)*,*(*) 

T*Q 

4//n-(z)n-(z)ff(z)^|^(z)^(z) . 

T/ze variation with respect to the components Uj of Tl q can easily be 
computed to give 

SH , ■ N 1 • • oflP (z) 

'z) = X h (IF (z)) + —5 v Ik (z) 



<5IIj m dt 

so that the first set of equations holds. For the other set of Hamilton's 
equations, we first note that the second integral in Hq may be written, 
up to a divergence, as 

2 ffdlljiz) , ■ <9 2 iY(q,q) 

T*Q 

and the derivative of this with respect to IF (z) gives 
aj >, (q) aUj (z) a /, [ BO,, (£)„„.. ,dK(q,® 



c i^^r + a? e y ^r n (i » * (z> ] . (55) 



Similarly, we compute the derivative of the last term in H with respect 
to the components IP (z) of U p to be 



dqidq 1 dqi \ 2 



T*Q 

Collecting these results we find 

SHn d 2 (bAa) . d<& (a) 

_ <m j (z) g$(q) 
dt dqi 
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where $ (q) is the function consisting of two gradient terms in Eqs. ( T55j) 
and (ESJ). In passing to the density formulation the additional gradient 
term is acted on by ^-divergence which identically vanishes because it 
is a function of q only. Thus, the Vlasov equations resulting from the 
canonical flow of the quadratic Hamiltonian in Eq. / f53j) and the momen- 
tum Vlasov equations are the same. ■ 

The canonical Hamiltonian system in proposition (!T4|) can be re- 
garded as a manifestation of the fact that, the addition of an exact 
one-form to 11^ does not affect the Vlasov equation in /. In fact, there 
are infinitely many such flows in momentum formulation due to symme- 
tries in the definition of momentum variables. To this end, we observe 
the following symmetries associated with the momentum formulations. 
The quadratic Hamiltonian Ho(Uid) is invariant under shifts in velocity 
together with a reflection in position 

n, ^ Ui + 2mX h {Yi l ) , U l ^ -rf. (57) 

Both the Hamiltonian and the canonical symplectic structure are invari- 
ant under arbitrary diffeomorphisms in the variables n\ The definition 
f fT6|) of the plasma density / is also invariant under 

n 4 ^ Ui + (V p x A p ), , IT ^ IT + (V, x A q y (58) 

for arbitrary vector functions A 9 (z) and A p (z). 

Regarding the components of U q as momentum variables in the 
canonical structure of proposition (114j) . Eqs. (140 j) becomes inverse Legen- 
dre transformations to be solved for the momenta. Then, the Lagrangian 
functional 

involving the velocity dYI p /dt shifted by the term —X h (II p ), gives the 
Euler-Lagrange equations 

rr(z) + 2X h (ti i ( z )) + x 2 h (u\z)) + -s ij ^^-u k (z) = o 

which can also be obtained from Eqs. ( 13 9 j) and (|4*0|) by eliminating the 
variables U^. 

7 Comparison with 2D Incompressible Fluid 

The motion of an incompressible fluid in a two dimensional region Ai C 
1R 2 is described by volume (area) preserving diffeomorphisms. The gen- 
erators are divergence-free vectors. Since the dimension is two, these are 
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equivalent to canonical diffeomorphisms and the generators are canon- 
ical Hamiltonian vectors of the form v = — where the Hamiltonian 
function ip is the stream function. The curl of v or, equivalently, the 
Laplacian of the stream function ip is the vorticity. Geometrically, it is 
a two-form given by 



d (v ■ dz) = V 2 z tpdq A dp, z = (q, p) (59) 



or, as a function u = V z ■ Q M (v ■ dz). 

The dynamics is governed by the Euler equation in vorticity form 

du 



dt 



where { , } M is the canonical Poisson bracket on Ai. This is the Lie- 
Poisson equation on the dual of the Lie algebra of the group of volume 
preserving diffeomorphisms for the Lie-Poisson structure 

W«).*H>^ = /„«<■) {^.^} A (•») 



with the Hamiltonian functional 



H= l -J v 2 d 2 z= l -J (V^) 2 d 2 z = ~Jl*> d 2 z (61) 

where a divergence term in the last expression is omitted. Alternatively, 
defining the Green's function solution 

i> (z) = - J K(z\z')u (z') d 2 z' (62) 

for the equation to = V z ip we have 

H = ^J J oj (z) K(z\z')u (z') d 2 zd 2 z'. 



See references [16],[28j from which we extract the above summary, and 
[31]- [35] for more on fluid motions. 

We observe that the quadratic Hamiltonian functional for incom- 
pressible fluid is a direct consequence of the definition of the dual of the 
Lie algebra by a metric. In this case, the (weak) non-degeneracy of the 
pairing is the same as the non-degeneracy of the metric and the Lie al- 
gebra can be identified with its metric dual. On the other hand, the Lie 
algebra g of Hamiltonian vector fields and its dual q* are L 2 — orthogonal 
in TT*Q. 
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The moment um-Vlasov equations in components of expresses the 
evolution of a volume cell, that is the density /, in the phase space T*Q 
in terms of its boundaries, that is, surfaces of the momenta 11^. This 
interpretation was first given by Ye and Morrison in [36J for the Clebsch 
variables (a, 0) defined by {a,/3} T * 2 = /. In the present context, they 
form a non-closed one-form ad (3 and can be identified with H^. 
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Inspired from the relation u = between the vorticity and Hamil- 
tonian functions of 2D fluid, we can establish a similar relation be- 
tween the plasma density function / and the Hamiltonian function h 
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of particle motion. The Hessian of h can be considered to be a map 
Hess(h) : TT*Q — > T*T*Q which is non-degenerate if the potential 
function 0j is non-degenerate. Let X e TT*Q, and define the vector 
field 

Y = tt i TtQ oHess(h)oX . (63) 

If X is a Hamiltonian vector field with a Hamiltonian function which is 
at least quadratic in momenta, then Y is not Hamiltonian. In particular, 
if we choose the Hamiltonian function to be h and identify Y with nj| d , 
then we get the Poisson equation. In other words, the Poisson equation 
in plasma resembles the relation between vorticity and stream functions 
of 2D fluid. This relation may also be described by assuming a non- 
degenerate Lagrangian functional on TT*Q which is yet to be found. 
Suppose we have a Lagrangian I on TT*Q quadratic in the velocities 
(q, p). We can introduce the momenta which reads 

U * = ITi = -P* > n = — = -e5 J —j 

oq l m dpi oqi 

for the special choice h of the Hamiltonian function. Then, the definition 
of plasma density in terms of momentum variables gives / = 1/m + 
eVg(j)f. Note also that with a rescaling of m and a redefinition of / we 
can write / = tr(Hess(h)). 

8 Conclusions 

Gauge symmetries of the Hamiltonian motion of the plasma particles 
leads to the kinematical constraint described by the Poisson equation. 
Thus, the Poisson type equations naturally arise in kinetic theories of 
particles moving in accordance with a canonical Hamiltonian formula- 
tion. Moreover, this implies that the true configuration space appro- 
priate for the dynamical formulation of the collisionless plasma mo- 
tion in Eulerian variables is the semi-direct product space J-(Q) (D 
Dif fcan(T*Q). This is well suited for a geometric understanding of 
the limit c i— > oo of the Maxwell- Vlasov equations. 

The formulation of dynamics in density variable is obtained by fur- 
ther reduction of momentum- Vlasov equations by the symmetry defining 
the gauge equivalence classes of momentum variables. The gauge alge- 
bra is, as a vector space, shown to be the same as $j but with an action 
different from the coadjoint action. The Eulerian velocity and momenta 
are complementary in the vector space TT*Q. Obviously, this and other 
geometric properties disappear upon identification of q and g* with func- 
tion spaces J^T^Q) and Den(T*Q), respectively. As an example, the 
function hf and the density / appear symmetrically in the Hamiltonian 
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functional of the Lie-Poisson structure whereas the corresponding vari- 
ables Xh f and Hid are complementary in the sense that f^» Q ({j) and 
0* decompose the space of one-forms on T* Q into spaces of exact and 
non-closed one-forms, respectively. The momentum formulation clarifies 
the geometric relation between the motions of plasma particles and the 
Lie-Poisson description of dynamics (23] ■ We expect the space TT* Q be 
important for Euler-Poincare formulation of dynamics [7] . [37] . [38] . and 
for application of Tulczyjew construction for Legendre transformation 
[20] from Lie-Poisson formulation. 
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